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ABSTRACT 

Accurate determination of the rotation rate in the radiative zone of the sun from helioseismic 
observations requires rotational frequency splittings of exceptional quality as well as reliable 
inversion techniques. We present here inferences based on mode parameters calculated from 2088- 
days long MDI, GONG and GOLF time series that were fitted to estimate very low frequency 
rotational splittings {y < 1.7 mHz). These low frequency modes provide data of exceptional 
quality, since the width of the mode peaks is much smaller than the rotational splitting and 
hence it is much easier to separate the rotational splittings from the effects caused by the finite 
lifetime and the stochastic excitation of the modes. 

We also have implemented a new inversion methodology that allows us to infer the rotation 
rate of the radiative interior from mode sets that span I = 1 to 25. Our results are compatible 
with the sun rotating like a rigid solid in most of the radiative zone and slowing down in the 
core {t/Rq < 0.2). A resolution analysis of the inversion was carried out for the solar rotation 
inverse problem. This analysis effectively establishes a direct relationship between the mode set 
included in the inversion and the sensitivity and information content of the resulting inferences. 
We show that such an approach allows us to determine the effect of adding low frequency and 
low degree p-modes, high frequency and low degree p-modes, as well as some g-modes on the 
derived rotation rate in the solar radiative zone, and in particular the solar core. We conclude 
that the level of uncertainties that is needed to infer the dynamical conditions in the core when 
only p-modes are included is unlikely to be reached in the near future, and hence sustained efforts 
are needed towards the detection and characterization of g-modes. 

Subject headings: solar interior rotation, inversions 
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Introduction 



solar radiative interior, as well as the level of accu- 
racy that helioseismic data should have to resolve 
the solar core. 



Over the past decade, increasingly accurate 
helioseismic observations from ground-based and 
space-based instruments have given us a reason- 
ably good d escription of the dynamics of the solar 

inter ior (e.g. lSchou et al,lll998tlThompson. ChristensendTa lsgaUidslMtHiA.P&iilEoaiiaifclinear rotational helio 



2. Theoretical background 



Helioseismic inferences have confirmed that 
the differential rotation observed at the surface 
persists throughout the convection zone. There 
appears to be very little, if any, variation of the 
rotation rate with latitude in the outer radiative 
zone (0.4 > t/Rq > 0.7). In that region the ro- 
tation rate is almost constant (« 430 nHz), while 
at the base of the convection zone, a shear layer 
— known as the tachocline — separates the region 
of differential rotation throughout the convection 
zone from the one with rigid rotation in the radia- 
tive zone. 

Despite the large scatter amongst the rotational 
splittings th at are sensitive to the solar core (see 
discussion in Eff-Darwi ch. Korzennik. fe Jimenez-Revesi 



2002), we can rule out an inward increase or de- 



crease of the solar internal rotation rate down 
to r/RQ w 0.25, by more t han 20% of the sur- 
face rate at mid-latitude (IChaol in et all 11999: 



Eff-D arwich. Korzennik. fc Jimenez-Revesi 



2002: 



Couvidat et aLl2*0*03 ). This is in clear disagree 



seismic inversion methodologies is the functional 
form of the perturbation in frequency, Au n i m , in- 
duced by the rotatio n of the sun, tl(r, 9), and given 
by (see derivation in lHansen et a/.lll977l ): 

Av nlm = — [ [ K n i m (r,6)tt(r,0)drdO±e n im 
2tt Jo Jo 

(1) 

The perturbation in frequency, Av n £ m , with the 
observational error, e n im, that corresponds to the 
rotational component of the frequency splittings, 
is given by the integral of the product of a sen- 
sitivity function, or kernel, K n £ m (r,9), with the 
rotation rate, Q(r,9), over the radius, r, and the 
co-latitude, 9. The kernels, K n g m (r, 9), are known 
functions of the solar model. 

Equation [1] defines a classical inverse problem 
for the sun's rotation. The inversion of this set 
of M integral equations - one for each measured 
Avnim - allows us to infer the rotation rate profile 
as a function of radius and latitude from a set of 
observed rotational frequency splittings (hereafter 
referred as splittings). 

The inversion method we use is based on the 



ment with the theoretical hydrodynamical models 
that predict a much faster rotation in the solar 
core, na mely 10 to 50 times faster than the surface 
rate ( e.g. lThompson. Christensen-Dalsgaard. Miesch. fc Tjkgala rized Least-Squares methodology (RLS). 

The RLS method requires the discretization of 
the integral relation to be inverted. In our case, 
Eq. [1] is transformed into a matrix relation 



2003). 



Mo re recently, iGarcia et all ( 2004f ) and Korzennikl 
J2005J) have independently developed new mode 
fitting procedures to improve the quality and pre- 
cision of the characterization of the modes that 
are sensitive to the rotation in the solar core. 
By using very long time series — spanning nearly 
six years of observati ons — coll ected with the MDI 
(|Scherrer et aZj|l995h. GONG (lHarvev et aZj|l996h 
and GOLF (jGabriel et aLlll995l ) instruments they 
have measured rotational splittings for modes with 
frequencies as low as 1.1 mHz. 

We present here an attempt to constraint the 
radial and latitudinal distribution of the rotation 
rate in the radiative interior through the inversion 
of a combined MDI, GONG & GOLF data set. 
We also attempt to establish the sensitivity of he- 
lioseismic data sets to the dynamics of the inner 



D = Ax 



(2) 



where D is the data vector, with elements Av n £ m 
and dimension M, x is the solution vector to be 
determined at N model grid points, A is the ma- 
trix with the kernels, of dimension M x N and e 
is the vector containing the corresponding obser- 
vational uncertainties. 

The RLS solution is the one that minimizes 
the quadratic difference x 2 — \Ax — D\ 2 , with 
a constraint given by a smoothing matrix, H = 
G T G, introduced in order to lift the singular na- 
ture of the problem (see, for additio nal details, 
Eff-Darwich fc Perez-Hernandez 19971 ). The ma- 



trix G represents the first-order discrete differ- 
ential operator, although it will be shown below 
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that the inversion technique we have developed is 
to a first order approximation independent of the 
choice of G. The general relation to be minimized 
is 

S(x) = (Ax - D) T {Ax - D) +-/xHx (3) 

where 7 is a scalar introduced to give a suitable 
weight to the constraint matrix H on the solution. 
Hence, the function x is approximated by 



Xcst 



(A T A + 1 H)- 1 A T D 



(4) 



Replacing D from Eq. [2] we obtain 



^cst 



(A T A + 1 H)- 1 A T Ax^ t Rx (5) 



hence 



R = (A r A + 1 H)- l A T A 



(6) 



The matrix R, that combines forward and in- 
verse mapping, is referred to a s the resolution or 
sensitivity matrix (Fricdel 2003). Ideally, R would 
be the identity matrix, which corresponds to per- 
fect resolution. However, if we try to find an in- 
verse with a resolution matrix R close to the iden- 
tity, the solution is generally dominated by the 
noise magnification. The individual columns of R 
display how anomalies in the corresponding model 
are imaged by the combined effect of measurement 
and inversion. In this sense, each element Rij re- 
veals how much of the anomaly in the j th inver- 
sion model grid point is transferred into the i th 
grid point. Consequently, the diagonal elements 
Ra states how much of the information is saved in 
the model estimate and may be interpreted as the 
resolvability or sensitivity of Xi. We defined the 
sensitivity \ of the grid point X4 to the inversion 
process as follows: 



Ri. 



Yfj=l R K 



(7) 



With this definition, a lower value of A^ means a 
lower sensitivity of Xi to the inversion of the solar 
rotation. We define a smoothing vector W with 
elements w; = A^ 1 that is introduced in Eq. 0]to 
complement the smoothing parameter 7, namely 



(A A + jWH)~ A D 



(8) 



Such substitution allows to apply different regu- 
larizations to different model grid points Xi whose 



sensitivities depend on the data set that are used 
in the inversions. In this sense, the inversion is a 
two steps process: first R is obtained from Eq. [5] 
for a small value of the regularization parame- 
ter 7. Then, the smoothing vector W is calcu- 
lated through Eq. [7] and the inversion estimates 
are obtained through Eq. [8] A set of results can 
be calculated for different values of 7, the opti- 
mal solution being the one with the best trade- 
off between error propagation and the quadratic 
difference \ 2 = lAzw, — -PI 2 as introduced in 



Eff-Darwich fc Perez-Hernanded (| 19 9 71 ) 



In this paper we show how to use the matrix 
R to study the sensitivity of helioseismic data sets 
to the rotation rate of the solar interior. Conse- 
quently, we present a theoretical analysis of the 
effect of adding low frequency and low degree p- 
modes, high frequency and low degree p-modes, 
and g-modes on the rotation rate of the solar core 
derived through numerical helioseismic inversion 
techniques. 

3. Observational mode parameters and 
inversion results from 2088-days long 
MDI, GOLF and GONG time series 

The work presented here is based on rota- 
tional frequency splittings measured from obser- 
vations by the GONG ground-based network and 
the MDI and GOLF experiments on board the 
SOHO spacecraft. All rotational splittings were 
computed from 2088-days long time series, start- 
ing April 30th 1996 and ending January 17th 2002, 
as summarized in Table [U 

The three data sets KM, KG & GG (see Tabled] 
for explanations) contain for the first time very 
low frequency rotational splittings (y < 1.7 mHz). 
These low frequency modes provide data of excep- 
tional quality, since the width of the mode peaks 
is much smaller than the rotational splitting. It 
is therefore much easier to separate the rotational 
splittings from the effects caused by the finite life- 
time and the stochastic excitation of the modes. 
The data set SM (see again Table [1} was obtained 
by averaging all the data sets result ing from fit- 
ting t he 72-days long MDI time series (jSchou et al 



1998) that overlap the April 30th 1996 to January 
17th 2002 period. The averaging process reduces 
significantly the number of I < 8 modes in that 
data set. 
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Since these data sets have been calculated from 
different time-series and peak-fitting techniques, 
one can expect some differences among the dif- 
ferent data sets. When using different time-series, 
the mode parameters can be affected by the chang- 
ing solar activity cycle. Moreover, fitting tech- 
niques can give different results if they are applied 
to eit her individual peaks or to ridges ( Korzennikl 
2005). Differences between MDI, GOLF and 



GONG can also arise from systematics introduced 
by the merging process used by GONG to ob- 
tain single time-series from multiple stations lo- 
cated worldwide. Differences may also come from 
the different spatial filters and leakage matrices 
that ar e used to isolate t he signal of an indiv idual 
mode (jKorzennikl 120051 IChaplin et al\\200$ ). In 
any case, we combined the various data sets in a 
single set - following the prescription described in 
Table [2] Our newly developed inversion method- 
ology was applied to the combined set to infer the 
rotation rate in the solar interior. 

Figure 1 shows the observational frequency 
splitting uncertainties of the combined data set 
as a function of radial order and degree, whereas 
Figs. 2 and 3 show the splittings uncertainties for 
sectoral modes as a function of a proxy of the inner 
turning point of the modes, i/v, and as a function 
of frequency, v, respectively. These plots clearly 
illustrate the well known and challenging fact that 
only a small number of modes penetrate the solar 
core and that the largest uncertainties are asso- 
ciated with these modes. Indeed, the combined 
data set does not include low degree high fre- 
quency modes (i.e., £ < 4 and v > 2.2 mHz), since 
at higher frequencies unwanted bias appears in the 
estimated splittings due to the difficulty in sep- 
arating the effect of rotation al splitting from the 



limited lifetime of the mo des (jAppourchaux et al. 
\200d IChaplin et aLlliooi l. 

The inversion of the combined MDI-GOLF- 
GONG data (see Fig. 4) confirms that the well- 
known differential rotation observed at the sur- 
face persists throughout the convection zone. Al- 
though only £ < 25 modes were used in the inver- 
sion, it was possible to infer the rotation rate in 
the convection zone as the result of the exceptional 
quality of the low frequency splitti ngs (y < 1.7 
mHz) obtained by iKorzenni 2 (|2005[) . The differ- 
ential rotation changes abruptly at approximately 
0.7 Rq to rigid rotation throughout the radiative 



zone. The radial distribution of the rotation is ap- 
proximately flat, at a rate of ss 430 nHz, decreas- 
ing below 0.2i?Q. The tendency below 0.15i?o 
is not real and results from extrapolation of the 
trend seen at larger radii, as explained in the fol- 
lowing section. 

4. Sensitivity analysis for the inversion of 
the solar radiative interior. 

A theoretical analysis was carried out in or- 
der to determine the effect of different low degree 
mode sets on the derivations of the solar rotation 
rate of the inner radiative interior. Four different 
artificial data sets, hereafter referred as Ai to A4 
were calculated using Eq. [1] and an artificial ro- 
tation rate OaC?") &) that is shown in Fig. 5. The 
different artificial data sets correspond to differ- 
ent mode sets and/or uncertainties, as explained 
in Table 3. The observational uncertainties (stan- 
dard errors) were taken from the combined mode 
set used in the previous section, whereas the noise 
added to the artificial data was calculated from 
normal distribution with the observed uncertain- 
ties. The Ai data set contains the same mode 
set than the combined MDI-GOLF-GONG data 
set. Errors for g modes were arbitrarily set to 6 
nHz (the mean of the observational uncertainties 
for the acoustic mode splittings), since at present 
there are not reliable estimates for the uncertain- 
ties of g-modes frequency splittings. In any case, 
we are interested in the behavior of the inversion 
methodology when g modes are added, rather than 
in the results of the inversion for different values 
of the splittings and the observational errors. 

The sensitivity vector, A, was computed for the 
four artificial data sets, as illustrated in Figs. 6 
and 7, where the sensitivities, Ai, for the rotation 
rate in the solar interior are presented as a func- 
tion of the radius, for each artificial data set at 
the equator or for several latitudes for a given set. 
The data sets Ai and A2 are significantly less sen- 
sitive to the rotation of the solar core than the 
other sets. Although the A2 set includes the same 
high frequency modes as set A3, the errors of the 
A2 set are significantly larger and hence the sen- 
sitivities do not differ from those obtained for the 
Ai set. The addition of two g-modes (in set A4) 
increases significantly the sensitivity to the solar 
core. However, it is important to notice that even 
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with the addition of g-modes, the sensitivity at the 
solar core varies with the latitude, as illustrated in 
Fig. 7. For all sets, the sensitivities to the equa- 
torial regions of the solar core are larger than the 
sensitivities at other latitudes. 

The effect on the sensitivity vector of the choice 
of the smoothing matrix H is presented in Fig. 8, 
where the equatorial sensitivities Ai for the first, 
second and third-order discrete differential oper- 
ators G are shown. The larger the order of the 
operator G the larger the sensitivity, except near 
the edges, and in particular the core. The smooth- 
ing vector W is obtained from the sensitivity vec- 
tor A and hence, the regularization constraint will 
not only depend on the mode set, but also on the 
shape of G. 

The choice of the number and spatial distribu- 
tion of the model grid points, N, is an important 
aspect of the inversion process. In non-adaptive 
regularization inversions, decreasing the number 
of grid points is in itself a form of regularization. 
The inversion procedure described here will also 
adjust to the distribution of grid points. This is 
illustrated in Fig. 9, where we show that the vari- 
ation of the inversion sensitivities (and hence the 
regularization constraints) is not constant with ra- 
dius and latitude when the number of grid points 
are changed. In this sense, the a priori choice of 
the distribution of model grid points will not con- 
strain the inversion results. 

The inversion methodology developed for the 
work presented here differs from standard RLS 
techniques by introducing a smoothing vector 
W. The purpose of this vector is to avoid over- 
smoothing the inversion solution in certain regions 
of the solar interior and hence loose valuable in- 
formation. This is illustrated in Fig. 10, where 
two different inversions of the same data set are 
presented, namely a standard RLS inversion (no 
vector W is added) dotted line, and the newly 
developed RLS inversion, solid line. The standard 
RLS inversion tends to over-smooth and hence 
to assign unrealistic low errors to the estimated 
rotation in the convection zone to get a stable 
solution in the radiative regions. However when 
W is added to the inversion procedure, the over- 
smoothing problem in the convection zone is mit- 
igated, since the sensitivities Xi in the convection 
zone are larger and hence, the smoothing coeffi- 
cients Wi are lower. The standard RLS technique 



would assign the same smoothing coefficients to 
all model grid points in the inversion. 

The correspondence between the sensitivity 
analysis shown in Figs. 6 and 7 and the infor- 
mation contained in the resolution matrix R is 
illustrated in Figs. 11 to 13, where we show the 
resolution vector corresponding to the estimate at 
the equator for r, = 0.06i? Q . T his corresponds to 
the av eraging kernel defined by iBackus fc Gilbert 
( 1970t ) for the Optimal Localized Averages tech- 
nique. In the ideal case, the resolution should be 
unity at the location where the solution is esti- 
mated and zero elsewhere. Only in the inversion 
of the A4 set is the largest amplitude of the res- 
olution vector centered at r, = O.O6i?0, and thus 
the result is reliable. The poor localization of the 
resolution at r, = O.O6i?0 for the inversions of the 
Ai to A3 set could not be improved by any inver- 
sion technique. However, this lack of resolution 
is taken into account by the sensitivity analysis, 
since it evaluates low sensitivities and hence as- 
signs large regularizations to those grid points. At 
larger radii both the location and the amplitude 
of the resolution are significantly increased, as il- 
lustrated in Fig. 14, where we show the resolution 
corresponding to the estimate at the equator for 
r, = 0.2i? Q obtained with the data set A3. 

The conclusions derived from Figs. 6 and 7 can 
also be drawn from the inversions of the data sets, 
as illustrated in Figs. 15 to 18. Figures 15 and 
16 show the inverted profile and error distribu- 
tion for the Ai set at several latitudes and demon- 
strates that there is good sensitivity to the latitu- 
dinal variation in the radiative rotation rate above 
r f» O.3i?0. Hence, the absence of differential ro- 
tation for the solar rotation rate in the radiative 
interior (Fig. 4) is real, not an artifact of the inver- 
sion procedure. The unrealistic flat rotation rate 
below r « O.15i?0 resulting from the inversions of 
the Ai and A2 sets is due to the lack of sensitivity 
of the mode set to that region. As a result the 
inversion extrapolates the trend of the solution at 
larger radii. There are not significant differences 
in the inversions of sets Ai and A2, although set 
A2 includes low degree and high frequency modes. 
However, the new information contained in the low 
degree and high frequency modes of set A2 is lost 
due to they large observational uncertainties. In 
all four cases, larger differences between the arti- 
ficial and the inverted rotation rates are seen at 
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higher latitudes, especially in the radiative inte- 
rior, as a result of the lack of sensitivity of the 
inversions to the polar regions. 

Only in the cases of sets A3 and A4 (see Fig. 17), 
was it possible to infer the main trends of the ro- 
tation rate below r w 0.15-R©. However, it was 
necessary to include, either, data with unrealis- 
tic small observational errors (set A 3 ), or, a cou- 
ple of (/-modes (set A4), modes that have yet to 
be unambiguously observed. The most likely way 
to reduce the observational uncertainties consists 
of increasing the length of the time series. Fig- 
ure 19 compares the observational errors for the 
I — 25 sectoral modes for five 728-days long data 
set s to the 2088-days long data set, all estimated 
bv iKorzennikl ((20oi . The formal observational 
uncertainties are proportional to the square root 
of the length of the time series, hence it would 
be necessary to observe for decades to reduce the 
observational uncertainties of the very low degree 
modes to the current levels of the I = 25 modes, 
all the while assuming that we can also reduce the 
residual bias in our current estimates of the low de- 
gr ee and high frequency splittings (see discussion s 



Appourchaux et al. 2000l ; Chaplin et al. 20061) 



5. Conclusions 

We have used for the first time a combined 
MDI-GOLF-GONG data set of rotational fre- 
quency splittings that covers the largest possi- 
ble frequency range, spanning from 1.1 mHz to 
3.9 mHz. This mode set was determined from 
2088-days long time series acquired by the MDI, 
GOLF and GONG instruments and ana lyzed inde- 
pende nt ly by several autho rs , namely Korzeniiikl 
(120051). iGareia et all (I2004I). iGellv et all (1 2002 k 
ISchou et all (|l998t ) and Ijimenez Revesl ([20001 ). 
Very low frequency splittings (y < 1.7 mHz) were 
included to improve the precision and the resolu- 
tion of the inversion in the solar interior. 

In order to optimally invert this unique data 
set, we implemented a new inversion methodology 
that combines the regularized least-squares tech- 
nique with the analysis of the sensitivity of the so- 
lution at all model grid point to the mode set being 
inverted. The inversion of the actual MDI-GOLF- 
GONG data set reveals that the sun rotates as 
a rigid solid in most of the radiative interior and 
slowing down below 0.2i? Q . 



The calculation of the sensitivity vector A of- 
fers a rapid and intuitive way of evaluating the 
sensitivity of helioseismic data to the dynamics 
of the solar interior, in particular in the core 
(r < O.25i?0). We conclude that with the present 
accuracy of the available splittings, it is not pos- 
sible to derive the dynamical conditions below 
r sa 0.2i? Q . This results from the relatively large 
observational uncertainties of the modes sensitive 
to the solar core, in particular the low degree and 
high frequency modes. The level of uncertainties 
that is needed to infer the dynamical conditions 
in the core when only including p-modes is un- 
likely to be reached in the near future, and hence 
sustained efforts are needed towards the detection 
and characterization of g modes. 
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Fig. 1. — Observational uncertainties of the fre- 
quency splittings, shown as a function of radial 
order n and spherical harmonic degree, £. Uncer- 
tainties obtained from 2088-days long time-series 
for the combined sectoral frequency splittings set 
are shown using a logarithmic gray scale. 



V 
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This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 



Fig. 2. — Observational uncertainties of sectoral 
frequency splittings as a function of the inner turn- 
ing point proxy ijv, for a selected set of degrees. 
The solid, dotted, dashed and dashed-dotted lines 
correspond to I = 25, 15, 5 and 1, respectively. 
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Table 1 : Description of the data sets used in this work. 



Instrument 


Degree 


Frequency (mHz) 


Reference 


Referred as 


MDI 


4 < I < 25 


1.1 < v < 3.9 


Korzennik 2005 


KM 


GONG 


4 < £ < 25 


1.1 < v < 3.9 


Korzennik 2005 


KG 


MDI 


5 < I < 25 


1.7 < v < 3.9 


Schou et al. 1998 


SM 


MDI 


1 < i < 10 


1.8 < v < 3.9 


Jimenez-Reyes et al. 2000 


JM 


GOLF 


1 < t < 3 


1.1 < v < 2.2 


Garcia et al. 2004 


GG 


GOLF 


e = i 


1.4 < v < 3.9 


Gelly et al. 2002 


BG 



Table 2: Description of the combined data set used, and a brief description of the combination procedure. 



Instrument 


Degree 


Frequency (mHz) 


Combination Procedure 


GOLF 


e = 1 


1.1 < v < 2.2 


Average of GG and BG 


GOLF 


^ = 2,3 


1.1 < v < 2.2 


Data from GG 


MDI 


1= 1,3 


2.2 < v < 3.9 


Data from JM 


MDI-GONG 


4 < I < 10 


1.1 < v < 1.8 


Average of KM and KG 


MDI-GONG 


4 < I < 10 


1.8 < v < 3.9 


Average of JM, KM, SM and KG 


MDI-GONG 


11 < I < 25 


1.1 < v < 1.8 


Average of KM and KG 


MDI-GONG 


11 < I < 25 


1.8 < v < 3.9 


Average of KM, KG and SM 



The MDI-GONG-GOLF set consists of the ensemble of these values, although the MDI data for £=1,3 and 
2.2 < v < 3.9 are only used to compute artificial data sets. 



Table 3: Description of the artificial data sets used to study the sensitivity of helioseismic data to the 
dynamics of the solar core. 



Frcq. range (mHz) 



Data set 


Note 


i=\ 


t = 


2,3 


i > 3 


Set Ai 




1 < v < 2.2 


1 < V 


< 2.2 


1 < v < 3.9 


Set A 2 




1 < v < 3.9 


1 < V 


< 3.9 


1 < v < 3.9 


Set A 3 


(1) 


1 < v < 3.9 


1 < V 


< 3.9 


1 < v < 3.9 


Set A 4 


(2) 


0.1 < v < 2.2 


0.1 < v < 2.2 


1 < v < 3.9 



^ Data set A3 differs from data set A2 in its uncertainties: the uncertainties for £ < 4 are set to be the ones for the 

sectoral splittings of £ — 25. 
^ One £ = 1 and one sectoral £ = 2 p-modes rotational splittings were added to the Ai data set. 
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Fig. 3. — Observational uncertainties of sectoral 
frequency splittings as a function of the mode 
frequency, v, for a selected set of degrees. The 
solid, dotted, dashed and dashed-dotted lines cor- 
respond to £ — 25, 15, 5 and 1, respectively. 



Fig. 4. — Inversion of the combined MDI-GOLF- 
GONG data set (see Table 2), plotted as a function 
of radius at several latitudes, 8. The solid, dot- 
ted, dashed and dotted-dashed lines correspond 
to the rotation rate for the equator, and at lati- 
tudes of 30°, 60° and 80°, respectively. Error bars 
are shown for the rotation rate at the equator, 
whereas the error distributions for other latitudes 
are shown in Fig. 12. 
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Fig. 5. — Artificial rotation rate for the solar in- 
terior shown as a function of radius, r, at differ- 
ent latitudes, 9. This artificial profile incorporates 
latitudinal variations in two zones to test the in- 
version sensitivity to latitude. The solid, dotted, 
dashed and dotted-dashed correspond to the rota- 
tion rate for the equator, and at latitudes of 30°, 
60° and 80°, respectively. 



Fig. 6. — Sensitivities of the solar rotation at 
the equator for sets Ai to A4, represented by 
solid lines, dotted lines, dashed lines, and dotted- 
dashed lines, respectively. 
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Fig. 7. — Sensitivities to set A4 of the solar ro- 
tation at different latitudes, namely the equator 
(solid line), 30° (dotted lines), 60° (dashed lines) 
and 80° (dotted-dashed lines). 



Fig. 8. — Sensitivities to set Ai of the solar rota- 
tion at the equator for three different smoothing 
matrices H, namely those corresponding to first 
differences (solid line), second differences (dotted 
lines) and third differences (dashed lines). 
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Fig. 10. — Inversions of set A3 as a function of 
radius at the equator. Results for the standard 
RLS inversion (no vector W) and for the newly 
developed inversion method are shown in dotted 
and solid lines, respectively. The artificial rotation 
rate used to calculate the input sets is shown as 
the thin solid line. 



Fig. 9. — Sensitivities to set Ai of the solar rota- 
tion at the equator for two different distributions 
of model grid points, namely 50 & 16 points (solid 
line) and 25 & 16 points in the radial and angular 
directions respectively (dotted line). 
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Fig. 11. — Plot of the resolution vector Rij, j = 
1, N, calculated for set Ai corresponding to the in- 
version estimated at the equator for j-j = 0.06i? Q 
(indicated by vertical dotted lines). The same re- 
sult was obtained for the inversion of A 2 . 
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Fig. 12. — Plot of the resolution vector Rij, j = 
1,N, calculated for set A 3 corresponding to the 
inversion estimated, as in Fig. 11, at the equator 
for r, = OMRr... 
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Fig. 13. — Plot of the resolution vector Rij, j = 
1,N, calculated for set A 4 corresponding to the 
inversion estimated , as in Figs. 11 & 12, at the 
equator for r, = 0.06i? Q . 
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Fig. 14. — Plot of the resolution vector Rij, j = 
1,N, calculated for set A 3 corresponding to the 
inversion estimated at the equator but for rj = 
0.20i?o (indicated by vertical dotted lines). 
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Fig. 15. — Inversion of set Ai as a function of 
radius at different latitudes, namely the equator 
(solid line), 30° (dotted lines), 60° (dashed lines) 
and 80° (dotted-dashed lines). The artificial rota- 
tion rate used to calculate the input set is shown 
as thin lines. 



Fig. 16. — Error distribution for the inversion of 
set Ai as a function of radius at different latitudes, 
namely the equator (solid line), 30° (dotted lines), 
60° (dashed lines) and 80° (dotted-dashed lines). 
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Fig. 18. — Error distribution for the inversions of 
set Ai to A4 plotted as a function of radius at 
the equator. Results for sets Ai to A4 are shown 
as solid, dotted, dashed and dotted-dashed lines 
respectively. 



Fig. 17. — Inversions of set Ai to A4 plotted as 
a function of radius at the equator. Results for 
sets Ai to A4 are shown as solid, dotted, dashed 
and dotted-dashed lines respectively. The artificial 
rotation rate used to calculate the input sets is 
shown as a thin solid line. 
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Fig. 19. — Observational uncertainties for = 25 
sectoral frequency splittings as a function of the 
mode frequency, v. The solid line represents 
the uncertainties for the 2088-days long data set, 
whereas the other line styles represent the uncer- 
tainties for the five 728-days long data sets. 
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